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The recently reported observation of VFTS 243 is the first example of a massive black-hole binary 
system with negligible binary interaction following black-hole formation. The black-hole mass (x10Mọ) 
and near-circular orbit (e ~ 0.02) of VFTS 243 suggest that the progenitor star experienced complete 
collapse, with energy-momentum being lost predominantly through neutrinos. VFTS 243 enables us to 
constrain the natal kick and neutrino-emission asymmetry during black-hole formation. At 68% confidence 
level, the natal kick velocity (mass decrement) is £10 km/s (<1.0Mo), with a full probability distribution 
that peaks when ~0.3M. were ejected, presumably in neutrinos, and the black hole experienced a natal 
kick of 4 km/s. The neutrino-emission asymmetry is $4%, with best fit values of ~0-0.2%. Such a small 
neutrino natal kick accompanying black-hole formation is in agreement with theoretical predictions. 
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Introduction.—Stars several times more massive than 
the Sun end their lives with the collapse of their iron cores. 
The explosion mechanism and the physics driving the 
formation of the compact object are still a matter of 
intense research [1-8]. In the delayed neutrino-driven 
mechanism, neutrinos revive the stalled shock wave, even- 
tually leading to a successful explosion. In this case, the 
stellar mantle is successfully ejected and the compact-object 
remnant is a neutron star (NS) in most cases. However, 
if the explosion mechanism fails, continuous accretion of 
matter onto the transiently stable proto-NS pushes the latter 
over its mass limit and a black hole (BH) forms. 

The extreme speeds of pulsars [9-11] indicate that NSs 
experience a natal kick during formation [12]. This natal 
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kick is attributed mainly to asymmetric mass ejection, but 
asymmetric neutrino emission may contribute [13-21]. 
Simulations of the collapse of massive stars show that 
hydrodynamical instabilities leading to large-scale asym- 
metries in the mass distribution—such as neutrino- 
driven convection [22-24], the standing accretion shock 
instability [25,26], or the lepton emission self-sustained 
asymmetry instability [27,28]—could also have impor- 
tant consequences on the natal kick of the compact 
remnant [16—-18,27]. In addition, fast rotation and the 
presence of magnetic fields at the moment of collapse are 
likely to affect the natal kick [16]. 

In the extreme scenario of complete collapse into a BH, the 
ejecta mass and natal kicks are thought to be very low 
(~1-10 km/s) [29-31]. In this case, mass-energy is lost 
via neutrinos and, to a lesser extent, gravitational waves [1,21]. 
This differs from the archetypical scenario in which aniso- 
tropic baryonic ejecta are the main carriers of momentum [32]. 

Observations of stellar-mass BHs in high-mass x-ray 
binaries (HMXBs) have been employed to constrain the 
impact of natal kicks of collapsing stars on the orbital 
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configuration of massive binary systems [33-35]. HMXBs 
are binary systems comprising OB-type stars that transfer 
matter via stellar winds [36] to their compact-object 
companion [37]. To date about 150 HMXBs have been 
discovered in the Milky Way, but only a handful of these 
are associated with BH companions [38,39]. 

Recently, spectroscopy [40—43] and astrometry [44,45] 
have allowed the first detection of inert (i.e., x-ray 
quiescent or dormant) stellar-mass BH binaries [46] in 
the field [47]. This population of BH binaries is especially 
intriguing. These objects have wider orbital periods than 
HMXBs and the stellar companions are well within their 
Roche lobes. This configuration, which implies little 
interaction following BH formation [49], makes inert 
BH binaries better probes of natal kicks than HMXBs. 
In this Letter, we present direct inference on neutrino natal 
kicks for the most massive inert BH detected to date: VFTS 
243 [41]. Our findings are summarized in Fig. 1. 

For BH binaries, the plausibility of the complete collapse 
scenario can be assessed because the impact on the orbital 
evolution of mass ejection (dM) and a natal kick (Vkick) 
during compact-object formation on binary star systems are 
well understood [67—69]. For example, near-instantaneous 
(quasi)spherically symmetric mass ejection during the 
stellar collapse results in a recoil (the Blaauw effect [70]) 
that leads to a systemic velocity of the binary as a whole, 
widens the orbit, makes it more eccentric, and can even 
disrupt the binary [71]. On the other hand, asymmetries 
in the ejected mass lead to a wider variety of configu- 
rations, potentially modifying the separation (increase or 
decrease), eccentricity, and inclination of the orbit. Such 
asymmetries could play an important role in the formation 
of BH-BH mergers [72], BH HMXBs [73], and inert BH 
binaries [74,75]. 

Properties of VFTS 243.—VFTS 243 belongs to the 
family of inert BH binaries recently discovered in the 
Milky Way [43] and the Large Magellanic Cloud [41] via 
spectroscopy (Fig. 2). It comprises a main-sequence O star 
with inferred mass of M, = 25.0+2.3M, and a BH 
companion with Mpy = 10.1 + 2.0Mọ, orbital period of 
P = 10.4031 + 0.0004 d, and eccentricity e = 0.017 + 
0.012, where the errors are the lo uncertainty inte- 
rvals [41]. The relatively high mass of the BH and the 
nearly circular orbit suggest that the system experienced 
complete collapse. With a stellar radius (R) well within 
the binary Roche lobe (fre = R/Rez 7% 0.33) [41], VFTS 
243 is a relatively wide BH binary (for example, in 
contrast, Cygnus-X1 is close to filling its Roche lobe 
with fry, Z 0.99 [76,77]). Moreover, the supersynchro- 
nously rotating massive star of VFTS 243 suggests that the 
effects of tides in synchronizing and circularizing the orbit 
are negligible [41]. 

Constraints on natal kicks.—In order to analyze the 
formation of VFTS 243, we use a semianalytic appro- 
ach [34] to calculate the probability that a circular 
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FIG. 1. Constraints on the natal kick from VFTS 243. This 
figure shows the probability distribution of matching the current 
orbital period, eccentricity, and systemic velocity of VFTS 243, 
given a certain mass loss (dM) and natal-kick magnitude (Vķick) 
following BH formation. In the heat map, darker (lighter) regions 
indicate a higher (lower) probability for a certain dM-v,;,. pair to 
reproduce the orbital configuration of VFTS 243. The blue 
shaded region, delimited by the solid blue vertical lines, shows 
the estimates for neutrino mass-loss decrements from stellar 
models of BH progenitors [66]. The side panels display the 
marginalized probability distribution functions (PDF, gray histo- 
grams) and cumulative distribution functions (CDF, red curves). 
With CDF(dM = 1.0M,) = CDF (vkick = 10 km/s) = 0.68, we 
conclude that solutions with low mass loss and small natal kicks 
are preferred. The marginalized distributions peak at around 
dM = 0.3M, and v,;.5 = 4 km/s, which is consistent with a 
solution where the mass decrement comes exclusively from 
neutrino emission. 


precollapse binary that received a natal kick during BH 
formation could reproduce the orbital configuration of the 
system as observed today. Models of the precollapse 
binary suggest it initially had a short (~1 d) orbital period 
and experienced mass transfer [41]. For such configura- 
tions, the tidal circularization timescale (~10* yr [78]) is 
significantly shorter than the time the stars spend on the 
main sequence (~10° yr [79]) prior to BH formation. 
Therefore, we assume the binary was circularized prior to 
collapse. Given that the current configuration is barely 
eccentric, it is reasonable to assume that the orbit was 
circular and the marginal eccentricity was induced during 
BH formation. Alternatively, the current eccentricity 
could be a small residual from incomplete circularization 
via tides prior to collapse. For the precollapse orbital 
configurations we assume independent, uniform wide 
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FIG. 2. Observed sample of field black-hole (BH) binaries with 
massive companions. For each binary, we show the BH mass 
(Mpp) on the abscissa, the stellar companion mass (M,) on the 
ordinate, and the eccentricity (e) through the color bar. We use 
circles to indicate high-mass x-ray binaries (HMXBs) and 
triangles for inert BH binaries (Inert BHB). The symbol of 
VFTS 243 is slightly larger for clarity. 


prior distributions on the orbital period between 
5 < P;/d < 15, on the natal-kick magnitudes between 0 < 
Vkick/ (km/s) < 200 and on the amount of mass ejected 
between 0 < dM/Mo, <5. We also assume that the direc- 
tion of the imparted natal kick is isotropic. We use the 
postcollapse orbital period, eccentricity, and systemic radial 
velocity as independent constraints [49], which we justify 
by the lack of correlation in the observationally inferred 
orbital parameters [41]. 

Figure | shows the results of our analysis. The heat map 
presents two distinct hot spots. In the limit where Vkick > 0, 
we have the solution for pure mass loss (Blaauw effect [70]), 
where the amount of ejected mass is linearly proportional 
to the newly acquired eccentricity of the binary. The 
other spot has dM > 0, i.e., a very small amount of mass 
is ejected asymmetrically in the frame of reference of the 
exploding star, and a natal kick makes the orbit eccentric. 
The marginalized one-dimensional probability distributions 
have modes at dM ~ 0.3Mo and vic, % 4 km/s. Explosions 
with moderate mass loss (21Mo) and natal kicks 
(210 km/s) can occur (see yellow in the heat map from 
Fig. 1), but require natal kick directions and magnitudes 
fine-tuned to balance the mass loss in order to form a low- 
eccentricity BH binary [74], as demonstrated by the diago- 
nal structure of Fig. 1. Under the assumption that the binary 
was circular prior to BH formation and has a small but 
nonzero eccentricity at present, there is no support for a 
solution with dM = Vķick = 0. 

In Fig. 1, the vertical blue band shows the estimated mass 
decrements associated with the total radiated neutrino 
energies of collapsing stripped stars [66]. Numerical sim- 
ulations of the collapse of massive stars without concomitant 


supernova explosion in 2D [31] and 3D [80] by the Garching 
group yield final neutrino-induced BH kick velocities of 
magnitudes similar to those favored by our analysis of VFTS 
243 (0—4 km/s), in particular also for a precollapse stellar 
model with similar mass. Contesting 3D simulations by the 
Princeton group found that the final neutrino-induced natal 
kick velocities of BHs formed via failed supernovae could be 
larger (7-8 km/s [81]), but still within our 1-o credible 
intervals. These quantitative differences may be connected 
to different life times of the NS prior to BH formation or to 
intrinsic differences of the neutrino asymmetries due to 
transport. 

Under the hypothesis of complete collapse, we perform 
a back-of-the-envelope estimate on long-term (œ 1 s) 
neutrino asymmetries during BH formation. We follow 
Refs. [14,29] and parametrize the linear momentum trans- 
ferred to the BH as a result of anisotropic neutrino emission 
as 


Mpuyvv kick = a, Efe = a,dM,c, (1) 


where q, is the neutrino emission asymmetry factor, V, kick 
is the neutrino natal kick, dM, is the mass decrement from 
neutrino emission, and c is the speed of light in vacuum. To 
estimate a, and v, ,ick We assume the following (see [49] 
for more details). We choose Mpy = 10Mo, similar to 
VFTS 243. We obtain the interval Et ~ 1-11 x 10°? erg 
from the stripped-star BH progenitors from Ref. [66]. 
Finally, the value of v,,icx is determined by obtaining 
the subset of the complete natal kick probability distribu- 
tion that is consistent with complete collapse. These 
assumptions result in an upper limit on the neutrino 
asymmetry of ~4%, where 0% implies spherically sym- 
metric neutrino emission (in the rest frame of the NS), and 
33% implies that the dipole component of the neutrino 
luminosity equals the monopole amplitude [80]. However, 
the asymmetries are between ~0-0.2% for the most likely 
solution of v, kiok & 3 km/s [49], though a solution with no 
natal kick is allowed. These values are in overall agreement 
with the estimates of the anisotropy parameter of the total 
neutrino emission in recent 3D core-collapse simulations 
with detailed neutrino transport, namely 0.45%-0.76% for 
successful supernova explosions [80] and 0.05%-0.15% in 
simulations of nonexploding models [80] (see, however, 
Ref. [21], which reports considerably higher values, prob- 
ably because of differences in the analysis). 
Discussion.—Our results provide evidence of neutrino- 
induced natal kicks and add strong support to the complete 
collapse formation scenario [41]. In general, the total natal 
kick of a compact object has a baryonic, neutrino, and 
gravitational-wave component. If matter ejection in a 
successful explosion takes place, it generally dominates 
the total kick for higher-mass progenitors [32,80,82,83], and 
gravitational waves contribute to the natal kick negligibly, 
since they carry orders of magnitude less energy with respect 
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to baryons and neutrinos (e.g., [31,83]). Our calculations 
provide constraints on the total natal kick and mass loss, 
which we find to be largely in agreement with mass loss 
exclusively through neutrino emission and an associated 
natal kick, rather than baryonic mass ejecta. Depending on 
the magnitude and direction, a natal kick from baryonic 
ejecta [75] could further reduce our limit on neutrino natal 
kicks. Alternatively, there could be a yet-to-be-determined 
physical mechanism that results in baryonic and neutrino 
kicks of a similar magnitude that are antialigned and almost 
completely cancel. 

Astrometric microlensing has been suggested as a method 
to constrain BH natal kicks [84,85]; as long as the observed 
velocity dispersion is much larger than the presumed 
neutrino kicks, the latter are unconstrained. The entire BH 
kick distribution cannot be completely explored with bound 
binaries [71], but systems similar to or more massive than 
VFTS 243 are ideal to explore and constrain low natal kicks. 
The broad implications of the natal kick that VFTS 243 
received during BH formation has been recently investigated 
in the literature [41], showing support for small [75] and 
moderate natal kicks [74]. The solution of very low natal 
kicks for VFTS 243 is in agreement with neutrino kick 
estimations from hydrodynamical simulations [31,86,87]. In 
particular, recent hydrodynamical simulations of core- 
collapse events yield a range of model-dependent neutrino 
kicks between ~0-4 km/s in 2D [31] and between 
0-3 km/s in 3D [80,87]. However, other simulations of 
BH-forming collapses produce neutrino natal kicks between 
30 km/s $ vy kick S 100 km/s [21], which would be parti- 
ally compatible with the tail of our distribution in Fig. | but 
may be difficult to reconcile with other existing theoretical 
constraints, if confirmed. 

For BH progenitors that have been stripped during a mass 
transfer episode, the complete collapse scenario suggests 
almost no baryonic ejecta and a mass decrement exclusively 
from neutrinos. A very massive proto-NS can lead to a 
total neutrino energy emission of E'*' ~ 1054 erg [66], which 
roughly corresponds to a neutrino mass decrement of 
dM? = 0.6Mo. Therefore, dM < 0.6Mo, and more real- 
istically 0.1 < dM/Mg, < 0.4 [49], is consistent with natal 
kicks induced from neutrino emission exclusively [66]. 
Alternatively, dM zZ 0.6M, implies at least some bar- 
yonic-mass ejecta. 

In addition to VFTS 243, several other inert BH binaries 
have been recently reported: HD 130298 [43], Gaia 
BHI [44,88], and Gaia BH2 [45], as well as inert BH 
binary candidates 2 MASS J05215658 + 4359220 [89,90], 
HD 96670 [40], VFTS 514 [42], and VFTS 779 [42], all of 
them presented in Fig. 2. Arguably, the most atypical 
systems in the sample are the astrometric BH binaries Gaia 
BH1 and Gaia BH2. With massive BHs (Mpgy = 10M), 
they are at wide orbital periods (= 100-1000 d), have low- 
mass (M,#1M,) companions and possibly formed 
dynamically [91,92], which could explain their large 


eccentricity (e ~ 0.5), putting them in a different class 
with respect to the other BH binaries, whose assembly is 
more in line with that of stellar binaries [46]. 

The sample of single line spectroscopic (SB1) inert BH 
binaries is diverse (Fig. 2). VFTS 243 is the most massive 
and least eccentric of SB1 BH binaries. HD 130298 is 
slightly less massive than VFTS 243 (with Mgy © 9Mọ and 
M,, ® 24M), has a similarly small Roche-lobe filling factor 
(frL = 9.26) and a similar orbital period (Pow ~ 15 d), 
but it is quite eccentric (e = 0.5). VFTS 514 is even less 
massive (with Mgy ¥5Mo and M, ~% 19M.) and is also 
quite eccentric (e ~ 0.4); with a long orbital period 
(Pow © 185 d), it is likely far from Roche-lobe overflow. 
VFTS 779 is the least massive of the SB! BH binaries (with 
Mpy %4Mo and M, ~ 14Mọ), has a near circular orbit 
(e x 0.02) and a long orbital period (x60 d). The spectral 
type of the star of VFTS 779 indicates it has evolved past the 
main sequence [42]; at this stage, it is developing a 
convective envelope that makes tidal circularization orders 
of magnitude more efficient than for main-sequence stars 
with radiative envelopes [49]. 

The sample of BH HMXBs [39,93] is broader in the 
mass parameter space, but otherwise is more homogeneous 
(Fig. 2). Most BH HMXBs have short orbital periods 
(Pow < 6d), small eccentricities (e < 0.1), and large 
Roche-filling factors (fp, Z 0.9). Large Roche-filling fac- 
tors in HMXBs are in agreement with the theory of wind- 
accreting x-ray binaries, which predicts a threshold of 
Fri Z 0.8-0.9 as a condition for observability [94]. Cygnus 
X-1 is the archetype of BH HMXBs, with Mpy x 21M, 
and M, x 41Mọ, a short-day orbital period (Po, % 6 d), a 
large Roche-lobe filling factor (fp, > 0.99), and a near- 
circular orbit (e ~ 0.02). While the eccentricities of Cygnus 
X-1, LMC X-1, and M33 X-1 are similar to that of VFTS 
243, both wind accretion onto the compact object [35,78] 
and tides could have played a role in circularizing these BH 
HMXBs. For massive, main-sequence stars with radiative 
envelopes, the dynamical tide is the dominant tidal dis- 
sipation mechanism [95]. Dynamical tides are not an 
efficient dissipation mechanism for VFTS 243 [42], the 
most compact inert BH binary detected to date. The 
circularization timescale via dynamical tides is a strong 
function of the separation [Tois « (a/R)?! [96], where a 
is the semimajor axis of the binary], which implies that 
tides are significantly less efficient for the SB1 sample than 
for the observed BH HMXBs. For example, doubling the 
orbital period of Cygnus X-1 would result in a similar 
period to VFTS 243 and would reduce the impact of tides 
by ~2 orders of magnitude. Moreover, HD 130298 has 
rather similar properties (masses, periods, and Roche- 
filling factors) to VFTS 243, but is rather eccentric; this 
demonstrates that dynamical tides are not efficient at 
circularizing the orbit at such large separations. 

Finally, we highlight that for BH binaries with M, > 
Mpy there seems to be a sharp drop in eccentricity for 
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systems with Mpy Z 10Mọ. We speculate that this drop 
could denote the transition between luminous, mass- 
shedding explosions and complete collapse [97—100]. The 
early theory of stellar collapse suggested that less massive 
cores eject some mass during BH formation, while complete 
collapse could only occur for cores of mass 711Mo [101]. 
However, a growing body of theoretical work shows 
evidence that the compact-object remnant mass function 
is a nonmonotonic function of mass [1,6,102—111] that 
could even be stochastic [112], dependent on initial rotation 
and composition [113], and that might be affected by binary 
interactions [114-116]. This means that more massive stars 
do not always result in more massive compact objects; some 
stellar cores above a certain threshold can result in mass- 
shedding explosions and form NSs [117]. Based on the sharp 
eccentricity transition between HD 130298 and VFTS 243 
(both SB1 BH binaries with similar masses and orbital 
periods) and the results presented in this Letter, we consider 
that helium cores with a mass of x10Mọ could undergo 
complete collapse. However, this does not rule out the 
possibility that complete collapse occurs for less massive 
stars, or that incomplete collapse returns at higher masses. 
Future observations should provide more precise values on 
the actual mass threshold for complete collapse. 

Conclusions.—We explore the effect of mass loss and 
natal kicks during BH formation in VFTS 243, the heaviest 
inert stellar-mass BH binary detected to date. Our results 
show that the marginalized distributions peak around 
dM = 0.3Mọ and vi, = 4 km/s, respectively; these val- 
ues are consistent with recent simulations of the stellar 
collapse predicting mass loss via neutrino emission exclu- 
sively [66], fully in agreement with the complete collapse 
scenario [102]. Under the hypothesis of complete collapse, 
we estimate the asymmetry in the neutrino emission to be 
~0%-0.2% and provide an upper limit of ~4%. We find 
that mass loss <1.0Mo and Vķick S 10 km/s are preferred 
at a 68% confidence level, in agreement with other studies 
in the literature [41,75]. Moreover, our results also accom- 
modate solutions with no natal kick and no asymmetries in 
neutrino emission. The progenitor of the BH component of 
VFTS 243 is likely a massive stripped star [42] that 
experiences negligible mass loss during BH formation. 
Following recent observations of inert BH binaries in the 
Local Group (Fig. 2), we suggest that complete collapse 
can occur for stars that end their lives with cores 
of 210M. 
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